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Abstract: We study thermodynamical properties of a fuzzy sphere in matrix quantum 
mechanics of the BFSS type including the Chern-Simons term. Various quantities are 
calculated to all orders in perturbation theory exploiting the one-loop saturation of the ef- 
fective action in the large- iV limit. The fuzzy sphere becomes unstable at sufficiently strong 
coupling, and the critical point is obtained explicitly as a function of the temperature. The 
whole phase diagram is investigated by Monte Carlo simulation. Above the critical point, 
we obtain perfect agreement with the all order results. In the region below the critical 
point, which is not accessible by perturbation theory, we observe the Hagedorn transition. 
In the high temperature limit our model is equivalent to a totally reduced model, and the 
relationship to previously known results is clarified. 
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1. Introduction 

Fuzzy sphere [1], which is a simple compact noncommutative manifolds, has been discussed 
extensively in the literature. One of the motivations comes from the general expectation 
that noncommutative geometry provides a crucial link to string theory and quantum grav- 
ity. Indeed Yang-Mills theories on noncommutative geometry appear in a certain low- 
energy limit of string theory [2]. There is also an independent observation that the space- 
time uncertainty relation, which is naturally realized by noncommutative geometry, can be 
derived from some general assumptions on the underlying theory of quantum gravity [3]. 
As another motivation, fuzzy manifolds may be used as a novel regularization method in 
quantum field theories [4]. 

In string theory, fuzzy spheres appear as D-branes in the presence of external fields [5]. 
In particular they appear as classical solutions 1 in the pp-wave matrix model [11], which is 
a generalization of the BFSS matrix theory [12] to the so-called pp-wave background [13]. 
Fundamental issues such as the stability of these solutions [10, 14, 15] and the spectrum of 
the fluctuations [14, 16-18] have been studied intensively. 2 

1 More general classical solutions such as a rotating fuzzy sphere are discussed in refs. [6-10]. 

2 Similarly, fuzzy spheres appear as classical solutions in matrix string theory [19] on a type IIA plane- 
wave background [20]. The spectrum around the fuzzy spheres is computed in ref. [21]. This theory is used 
to study the matrix big bang [22] . 
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Thermodynamical properties of the pp-wave matrix model have also been studied by 
various authors. In ref. [23] the free energy around the trivial vacuum, which corresponds 
to a transverse M5-brane [24] at zero temperature 3 , was evaluated at the one-loop level, and 
the Hagedorn transition was studied in detail. (See refs. [26,27] for a two-loop extension and 
ref. [28] for a review on this subject.) This calculation have been extended to more general 
vacua in refs. [29-31]. In all these works, however, a mass parameter, which parametrizes 
the deviation from the flat background, is assumed to be large so that higher loop effects 
can be neglected. 

In this paper we show that in fact it is possible to study the fuzzy sphere thermody- 
namics to all orders in perturbation theory. While the method can be applied to more 
general models including the pp-wave matrix model, here we demonstrate it in a simple 
model, which can be obtained by dimensionally reducing a 4d U(iV) gauge theory to Id. 
The model has been used recently to compute the mass gap in the theory of bosonic mem- 
branes [32]. The action contains the 3d Chern-Simons term representing the coupling to 
a background flux [5], which enables fuzzy spheres to appear as classical solutions. When 
the Chern-Simons coupling is switched off, it reduces to the "4d bosonic BFSS matrix the- 
ory" [33] . Studying matrix quantum mechanics at finite temperature 4 is itself an interesting 
subject [35,39-41], in particular, because of its relation to the black hole physics [42-44]. 

The method for the all order calculation has been developed in totally reduced models 
[45,46], which are motivated in the context of the type IIB matrix model [47]. The key 
observation is that, in the large- N limit, the effective action is saturated at one loop in the 
bosonic case [46], and at two loop in the supersymmetric case [45]. Through the Legendre 
transformation, one can obtain the free energy and various observables to all orders. (We 
emphasize that this is different from a typical situation in supersymmetric field theories, 
in which higher loop corrections simply vanish due to cancellation.) In the bosonic case, 
it has been confirmed that the all order results are in perfect agreement with the Monte 
Carlo results [48] obtained in the fuzzy sphere phase. These works have also been extended 
to four-dimensional fuzzy manifolds [49,50]. 

Similarly to the results in ref. [48], we find in the present finite-temperature system 
that a single fuzzy sphere becomes unstable at some critical a, the coefficient of the Chern- 
Simons term. This phenomenon occurs at any temperature, and we obtain explicitly the 
critical a as a function of the temperature. We also perform Monte Carlo simulation and 
confirm that the all order results for various observables agree very well with the Monte 
Carlo results above the critical a. In the region below the critical a, which is not accessible 
by perturbation theory, we observe the Hagedorn transition at some critical temperature. 
At high temperature our model is equivalent to a totally reduced model, which is analogous 
to the model studied in ref. [48]. We clarify the relationship to the results obtained there. 

The rest of this paper is organized as follows. In section 2 we define our model and 
discuss its classical solutions. In section 3 we show how one can perform the all order 

3 A fuzzy five-sphere solution was constructed [25] in a deformed plane-wave matrix model with an 
interaction term due to the 6-form potential. 

4 More generally, large- N gauge theory at finite temperature has been an active field of research [34-38] 
partly motivated from the gauge/gravity correspondence. 
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calculation in perturbation theory. In section 4 we compare the all order results with the 
Monte Carlo results. In section 5 we study the region in the phase diagram below the 
critical a, and show that the Hagedorn transition takes place. In section 6 we discuss the 
high temperature limit of the model. Section 7 is devoted to a summary and discussions. 

2. The model and its classical solutions 

The model we study in this paper is defined by the action 5 

S = N J^dttr^DtXiit)) 2 - (*),*,-(*)]) 2 + ^iae ijk X i (t)X j (t)X k (t)} , (2.1) 

where D t represents the covariant derivative D t = dt — i [A(t), ■ ] . The dynamical variables 
A(t) and X{(t) (i = 1,2,3) are N x N Hermitian matrices, which can be regarded as the 
gauge field and three adjoint scalars, respectively, in a Id gauge theory with the U(iV) 
gauge symmetry 

Xi(t) ^ g(t) Xi(t) g(rf ; A(t) - g(t) A(t) g(tf + i g(t) j f g(t? . (2.2) 

The Euclidean time t in (2.1) has a finite extent (3, which is related to the temperature T 
through (3 = 1/T, and all the fields obey periodic boundary conditions. The cubic term 
represents the Chern-Simons term, which is crucial for fuzzy spheres to become classical 
solutions. The a = case corresponds to the "4d bosonic BFSS model" studied in refs. [33]. 
The classical equations of motion can be obtained from the action (2.1) as 

(D t fX t = [X,, [Xj,Xi]] + i a e ijk [Xj,X k ] , (2.3) 
[X i ,D t X i ]=Q. (2.4) 

There are two types of static solutions. The first type is given by configurations with Xi(t) 
and A(t) being static and diagonal. The action vanishes identically for such configurations, 
and therefore all the diagonal elements are moduli parameters. The second type of solutions 
can be represented as 

s s 

X i (t) = 0(aL^ ) ®l fc7 ) , ^(t) = 0(l„ / ®i( / )) , (2.5) 
i=i i=i 

where represents the n-dimensional irreducible representation of the SU(2) algebra 
= ieij k L^ , and the parameters kj and nj satisfy ]C/=i n l • kj = N . The 
ki x fcj Hermitian matrices AS 1 ^ are arbitrary, and they represent the moduli parameters. 
For this type of classical solutions, the action is evaluated as 

7=1 

5 We could have replaced the overall factor of N in the action (2.1) by -p, where g represents the Yang- 
Mills coupling constant. Our choice would then correspond to setting the 't Hooft coupling A = g 2 N to 
unity. We do not lose any generality, however, since the model for arbitrary A can be readily obtained by 
rescaling X t -> \~ 1/z X h (3 A 1/3 /3, a A~ 1/3 a. 
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which becomes minimum for s = 1, k\ = 1, n\ = N. In this case the solution simply 
becomes 

X l (t) = aL ( i N) , A(t) = 0, (2.7) 

which represents a single fuzzy sphere with the radius p = \ aV N 2 — 1, since it satisfies 
Yl^=i(-^i) 2 = P 2 liv • ("Fuzzy" because of the non-trivial commutation relation among Xi.) 

Since the action evaluated for the fuzzy sphere type solutions (2.6) is proportional to 
a 4 , it is expected that the single fuzzy sphere (2.7), which gives the minimum action among 
those solutions, dominates the path integral at sufficiently large a. 



3. Perturbative calculation around the fuzzy sphere 
3.1 Exact effective action and the critical point 

In this subsection we calculate the one-loop effective action around a configuration Bi = 
kL^ N \ which reduces to the single fuzzy sphere solution for k = a. It is known that the 
effective action around a fuzzy sphere configuration is "one-loop exact" in the sense that 
higher order corrections vanish in the large- N limit [46,49,50]. From the effective action, 
we can obtain the critical coupling a c , below which the fuzzy sphere becomes unstable due 
to both quantum and thermal fluctuations. 

Let us first expand Xi(t) and A(t) around the rescaled single fuzzy sphere Bi as 

Xi(t) = Bi + Xi{t) , A{t) = + A{t), (3.1) 

where the fields Xi(t) and A(t) represent the fluctuation. Since the original action (2.1) 
has a gauge symmetry (2.2), we fix the gauge by adding the gauge-fixing term and the 
ghost term as 

Stotal = S + Sg.f. + S g h , (3.2) 
S e .t. = \n jdttr(d t A-i[Bi,Xi]) 2 , (3.3) 

S g h = N jdt\x(dtc- D t c- [Bi,c][Xi,c]) . (3.4) 
Plugging (3.1) into eq. (3.2), we obtain S tota ,i = S + Si + S 2 + S" 3 + S 4 , where 6 

So = - 1) Qk 4 - , (3.5) 

5 2 = N Jdttr{^Xi(-d? + K 2 £f)l, + \A{-Q 2 t + n 2 C 2 )A + c(-d 2 + k 2 £ 2 )c} , (3.6) 

5 3 = N Jdt tr( - [X^XjWB^Xj] + -iae ijk XiXjX k + c[B h [X t ,c]} 

-{[A,Bi}+id t Xi)[A,Xi]-id t c[A,c]) , (3.7) 

6 We have omitted a term — iV J dttx ([Bi, Bj] — iatijkBk) [Xi, Xj] in eq. (3.6), which does not contribute 
to the effective action at one loop. 



- 4 - 



while the linear term S\ and the quartic term S4 will not be needed in the following 
calculation. In eq. (3.6), we have introduced the adjoint operation £jM = [l[ N \m] on an 
N x N matrix M. Following the usual procedure, the effective action can be calculated 
as T(k) = r(°)(re) + T^(k), where the classical term is nothing but T(°\k) = So, and the 
one-loop term is given as 

T (1) (k) = lndet(-d 4 2 + k 2 C 2 ) (3.8) 

by performing the Gaussian integration over the fluctuation fields with the quadratic terms 
(3.6). When taking the determinant in eq. (3.8), we omit the zero mode corresponding to 
the constant mode proportional to the unit matrix. In order to diagonalize the operator 
(—d 2 +K 2 C 2 ) , we introduce the matrix analog of the spherical harmonics Y\ m (0 < I < N—l, 
—l<m<l), which obeys the orthonormal relations 

^{Y? m Y Vm ,)=5 lv 5 mm ,, Yl = {-l) m Y h _ m , (3.9) 



and has the following properties as a representation of the SU(2) algebra 

C-zYlm = mY tm , 

C 2 Y lm = 1(1 + l)Y lm , (3.10) 
£±Yim = T m)(l ± m + V)Y^ m±1 , 

where C± = C 1 ± iC 2 . Using the formula ]T£Li (} + S) = S ^F> we obtain 

N-i , „ s 

rW(«) = 2^(2; + l)lnisinh(^7Z(7Tl)) [ . (3.11) 
1=1 ^ ' 

Here we have omitted a ^-independent constant, 7 which is irrelevant for the following 
analysis. 

When we take the large- iV limit of the effective action, we have to scale the parameters 
a, (3 and k in such a way that the classical term T^°\k) and the one-loop term T^(k) 
become the same order. This motivates us to introduce the rescaled parameters 

a = N l ^a , $ = iV 2/3 /3 , k = N 1 ^k . (3.12) 



The sum over I in eq. (3.11) can be evaluated in the large- N limit with fixed (3 and k. Thus 
we obtain the exact effective action as 

1 .. 1 / 1 . 2 



Um o -^r(K) = ^(i« 4 -|a« 3 ) +^0k) = f(k;aj) . (3.13) 



The function $>(x) is defined as 

r-JV 



= \x-i ln(l - e») + 21n (sinh |) - -U 2 {e>) + -4jLi 3 (e») - -ic(3) , (3.14) 

O \ Zj / X X X 



7 This constant becomes relevant, e.g., when one compares free energy for different types of vacua [31]. 
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where the polylogarithm function Li n (z) and the Riemann zeta function C( n ) are defined, 
respectively, as Li n (z) = Y?k=i w and ((n) = 52kLi W- 

The local minimum of the effective action, which corresponds to the quantum fuzzy 
sphere, can be obtained by solving 

^f(k;aj) = (3.15) 

with respect to R in the region k ~ a. As we decrease a, we find that the local minimum 
disappears at some critical point a c , which depends on (3. The critical point d c obtained 
in this way is plotted against T = 1/$ in figure 2. In particular, the asymptotic behaviors 
of the critical point at the low T and high T limits are given by 



a c = < 



' 9 1 / 3 ~ 2.08 at f = , 

(UgLf) V4 ^2.48T V4 at T»l. (3 " 16) 



3.2 One-loop calculation of observables 

In this subsection we calculate the expectation values of the operators 



M = T^Ta I dteijktriXiXjXk) , 
o iv p Jo 

around the single fuzzy sphere (2.7) at one loop. Unlike the effective action, the expectation 
values do have higher-loop corrections, which shall be obtained in a resummed form in the 
next subsection. 

Let us decompose the fields into the background and fluctuations as in eq. (3.1), where 
we set k = a in this subsection. The expectation value (R 2 ) can be represented as 

{r2) = W I dttrB * + W J dttlB ^^ + w]3 fa**?®) • ( 3 - 18 ) 

The first term can be easily evaluated as 

^ fdttrBf = \a\N 2 -l) . (3.19) 

The second term can be evaluated at one loop using the cubic terms (3.7) as 

A JdttrBiiXiit)) = ^(/dttr(B i X i (t))/^tr([X J -(O,X fc (O][Si,^fc(*')])) 

-^(/*tr(5 i X i (t))/(^tr(c(*')[B J -jX,-(0,c(t')]])) o 



+^(/dttr( J B J A > J (t))Ktr(^(i / ),^]^(i , ),^(i / )]))„ ,(3-20) 



o 
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where the symbol ( • }o represents the expectation value using the quadratic terms (3.6) 
only. Eq. (3.20) can be evaluated by using the Wick theorem. The propagators can be 
derived from the quadratic terms (3.6) as 

((X i (t)) p9 (X j (t / ))rs)0 = SijApqrsit ~ t') , (3.21) 

((A(t)) pq (A(t')) rs )o = A pqrs (t - t') , (3.22) 
((c{t)) pq (c{t')) rs ) = A pqrs (t - t') , (3.23) 

where the indices p, q, r, s run over 1, • • • ,N and A pqrs (t — t') is defined as 

I 00 N ~ 1 f 1 (_i^m e 2nin(t-t')//3 
\ q rs(t-t') = j^ EE'E (2vrn//3) 2 + ^^ + l) (n " m)p<?(yZm)rS ' ^ 

n=-oo /=0 m=-l y /f J V ' 

The symbol Yl' implies that the zero mode is omitted by excluding I = for n = 0. Using 
the formula XmLi x ' 2 +n' 2 = ~ lx + lx c °th(x7r), eq. (3.20) can be evaluated as 

N—l 

A JdttrBiiXi®) = ^(2Z + l)^(ITl)coth(^V^+iy) • (3-25) 

The sum over / can be evaluated at large N as in (3.14) for fixed a and $, and it turns out 
that (3.25) is given by — ^$'((3 a). Since the third term of eq.(3.18) is suppressed at large 
N, (R 2 ) is obtained at one loop as 

lim ^(R 2 )^ = \o? - \&{fla) . (3.26) 

The expectation values of M and F 2 can be calculated in a similar way, but it is much 
easier to obtain them by making use of the fact that these operators appear in the action 
(2.1). The expectation values can therefore be rewritten as 8 

< M > = ik^ w{a ' m • (3 ' 27) 

Here the free energy W(a, (3) is defined by 

W(a,f3) = -In (^J[dX][dA]e- s ^j , (3.29) 

and at one loop it can be obtained from the effective action by simply replacing k by a. 
In the large- ./V limit with fixed a and (3, we get 

lim -L Wl „ loop (a,P) = -^/3a 4 + <f>0&) . (3.30) 

iv— >oo iv Z4 



8 Eq. (3.28) can be derived by introducing a source term in the action, and by absorbing it by rescaling 
the variables as ! h ^ 1,>3 t, Xi 1— » [iT^^Xi, A 1— > [i x l A A with an appropriate /i. Since the integration 
measure and the kinetic term in the action are invariant under this transformation, the free energy for the 
action with the source term can be obtained by simply rescaling a and (3. 
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Plugging this into (3.27) and (3.28), we obtain 



lim -j-< M >i-ioop = -Ia 3 + $'(/3d) , 



-4(^ 2 )i-ioop = \a A -2a&(pa) 

AT->oo 2 



lim 



(3.31) 
(3.32) 



3.3 All order calculation of observables 



In this subsection we exploit the fact that the effective action is saturated at one loop in 
the large- TV limit, and calculate the expectation values of the operators R 2 , M and F 2 to 
all orders in perturbation theory. The crucial point here is that the free energy and the 
effective action are related to each other by the Legendre transformation. Therefore, we 
can obtain the free energy by evaluating the effective action at its local minimum. Since the 
expectation values can be obtained by differentiating the free energy (for an action including 
an additional source term if the operator does not exist in the original action) , we can obtain 
the all order results for the expectation values in the large- N limit. This amounts to [46] 
keeping only the terms in the one-loop result that come from 1PI diagrams, and replacing 
a by the solution to eq. (3.15), which we denote as ko in what follows. Since the one-loop 
contributions to (M) and (R 2 ) come only from 1PR diagrams, the corresponding all order 
results are readily obtained from the classical results by replacing a by ko as 



1 



(^)all-ordcr = ^i^oY 



m 

1 1 3 

-77 ( M ) all-order = -~(«o) 



(3.33) 
(3.34) 



N x ' 6 

Let us next consider (F 2 ). Since the one-loop contribution to (F 2 ) includes both 1PI 
diagrams and 1PR diagrams, it is easier to obtain the all order result by using the relation 
(3.28). As explained above, the free energy is given to all orders in perturbation theory as 

1 



J™ ^all-ordcr^,/?) 
N—>oc iV z 



f(k ;a,P) . 



(3.35) 



When we differentiate Wall-order P) with respect to a and (3, we have to take into account 
that ko depends on a and 0. Thus we obtain 



lim ¥ {F )all-order 

N~*oo jY3 

5_ d l~,d 
—a— + -0—= 
6 da 3 dp 



4 

P 
4 

P 



f{k ;a,P) 



- a (— A— 

6 \da dk 



where the coefficients A and B are given as 



A 



dk __ 


d 2 f(k;a,P) 




( d 2 f(k;aj) 


da 


dkda 


k=KO 


y d 2 k 


__ dk _ 


d 2 f(k;a,P) 




( d 2 f(k;a,P) 


dp ~ 


dkdp 


K=KQ 


y d 2 k 



K=Kf)/ 



K=Kf)/ 



(3.36) 



(3.37) 



(3.38) 
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4. Comparison with Monte Carlo results 

In this section we compare the all order results obtained in the previous section with the 
results of Monte Carlo simulation taking the single fuzzy sphere (2.7) as the initial config- 
uration. The lattice formulation and the algorithm used for simulating the model (2.1) is 
the same as in ref. [51]. The lattice spacing a and the number of sites N t in the Euclidean 
time direction obey the relation N t a = (5. We have chosen these lattice parameters so that 
our results represent the continuum limit with sufficiently good accuracy. 9 

4.1 Boundary of the fuzzy sphere phase 
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Figure 1: The observables (i? 2 )/iV 4 / 3 , (F 2 )/N 2 / 3 , (M)/N and (|P|) arc plotted against a for 
T = 0.1. The dotted, dashed and solid lines represent the classical, one-loop and all order results, 
respectively. 



Let us first investigate how the observables (3.17) behave as a function of a. This, 
in particular, allows us to determine the critical a, below which the single fuzzy sphere 
becomes unstable. 

In figure 1 we plot the expectation values (R 2 )/N 4 / 3 , (F 2 )/N 2 / 3 , (M)/N against a 
for fixed T close to T = 0. Monte Carlo results show a discontinuity at a ~ 2.1, which 



9 More precisely, the lattice parameters are chosen to satisfy both a < e and N t > 10 at any temperature, 
where e = 0.02 is used for figures 1 and 3 (except for the right bottom panel), and e = 0.05 otherwise. See 
ref. [52] for an analysis on finite lattice spacing effects in a related model. 



-9- 



4.0 



o 



3.0 



2.0 




1.0 



0.0 



0.0 



1.0 



2.0 



3.0 



4.0 



5.0 



T 



Figure 2: The critical <S, below which the fuzzy sphere becomes unstable, is plotted against T. 
The open circles represent the results obtained by Monte Carlo simulation for TV = 16. The solid 
line represents the result obtained from the one-loop effective action in the large- N limit. 

agrees with the result (3.16) at T = 0. Above the critical point, Monte Carlo results for 
N = 16, 24 lie on top of each other as expected from perturbation theory, and they agree 
very well with the all order results given by (3.33), (3.34) and (3.36). We have also plotted 
the classical results and the one-loop results for comparison. It clearly demonstrate the 
existence of higher-loop corrections, which are included in the all order results. 
In Monte Carlo simulation we also calculate the Polyakov line 



where the symbol Pexp represents the path-ordered exponential. Results for (|P|) are 
shown in the right bottom panel of figure 1. We observe a gap at the same a as the other 
observables. The properties of the Polyakov line will be discussed later in more detail. 

From Monte Carlo simulations at various T, we obtain the critical a as a function of 
T, which is plotted in figure 2. We observe perfect agreement with the results obtained 
from the one- loop effective action in the large- iV limit. This confirms that the effective 
action is indeed saturated at one loop. 

We call the region above the critical line the fuzzy sphere phase, and the region below 
the critical line the Yang-Mills phase, following the terminology used in ref. [48]. The 
phase transition between the fuzzy sphere phase and the Yang-Mills phase continues to be 
of first order at any temperature, judging from the existence of discontinuity. In section 5 
we will see that the Yang-Mills phase is further divided into two phases by the Hagedorn 
transition. 

4.2 Temperature dependence of observables 

Next we investigate the temperature dependence of observables. In figure 3 we plot the 



expectation values (R 2 )/N 4 / 3 , (F 2 )/N 2 / 3 , (M)/N against f for a = 3.0. There is a gap 
at T ~ 2.0, as expected from figure 2. The all order results reproduce the T dependence of 




(4.1) 



the observables very well below the critical T. Thermal effects tend to shift the observables 
towards the values above the critical temperature. 
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Figure 3: The expectation values (R 2 ) /TV 4 / 3 , (F 2 ) /TV 2 / 3 , (M) /N arc plotted against f for a = 3.0 
and N = 16 . The dotted, dashed, solid lines represent the classical, one-loop, all order results, 
respectively. In the right bottom panel we plot (\P\) against T, focusing on the small T region, for 
a = 3.0 and N = 16 . The solid line represents a fit to eq. (4.2). 

In the right bottom panel of figure 3, we plot the Polyakov line (\P\) as a function of 
T. We have magnified the small T region in order to see how the Polyakov line decreases 
as T approaches 0. (Note that the scale of T in this plot is an order of magnitude smaller 
than in other plots in figure 3.) Our results can be nicely fitted to the behavior 



(|P|)=ex P (-4) , 



(4.2) 



which suggests that the system is in the "deconfined phase". The fitting parameter c = 
0.0063 corresponds to the energy increase caused by a single heavy "quark". From this 
figure we conclude that the center symmetry is always broken at T ^ 0. This statement 
needs some care, however. See footnote 11. 



5. Hagedorn transition in the Yang-Mills phase 

In this section we investigate the properties of the Yang-Mills phase. Perturbation theory 
is not applicable here, but Monte Carlo simulation continues to be a reliable method. 
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In figure 4 the Polyakov line 
is plotted against T. We find that it 
changes very rapidly at the temperature 
T ~ 1.1, which we denote as Th- Above 
Th, the data are clearly nonzero, and 
they have little dependence on N. Be- 
low Th, the data are consistent with 
decreasing as 1/N at large N. Thus our 
data suggest that the center symmetry is 
spontaneously broken at T > Th- This 
transition can be interpreted as the Hage- 
dorn transition [23,34], and the critical 
temperature Th is referred to as the Hage- 
dorn temperature in what follows. The 
value of Th is close to the result Th ~ 
A 1 / 3 obtained in ref. [33] 10 , where A is 
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Figure 4: The Polyakov line (\P\) is plotted against 
T for a = 0.0 and N = 16, 24, 32 . The dashed line 
represents the result obtained by eq. (6.7) in the 
high T limit. The solid lines represent the result 
including the next-leading order terms [57]. 



the 't Hooft coupling constant, which is set to unity in our analysis. (See footnote 5.) 



2.5 
2.0 



W A 1.5 

DC 



1.0 
0.5 



0.0 




N=16 
N=24 
N=32 
leading 
next-leading 



0.0 



0.5 



1.0 



1.5 
T 



2.0 



2.5 



3.0 




N=16 
N=24 
N=32 
leading 
next-leading 



1.5 

T 



3.0 



Figure 5: The observables (i? 2 ) and (F 2 ) are plotted against T for a = 0.0 and N = 16,24,32. 
The dashed lines represent the results obtained by eqs. (6.4) and (6.6) in the high T limit. The 
solid lines represent the results including the next-leading order terms [57]. 

In figure 5 we plot the observables {R 2 ) and (T 2 ) against T at a = for N = 16, 24, 32. 
The results for different N lie on top of each other, which implies a clear large- N scaling 
behavior. In the confined phase T < Th, we find that the results are independent of T. 
This can be considered as a consequence of the Eguchi-Kawai equivalence, 11 which states 

10 Note, however, that the lattice model studied in ref. [33] is written in terms of unitary matrices Ui(t) 
instead of Hermitian matrices Xi(t), and it agrees with our model only after replacing Ui(t) by exp(iaXi(t)) 
and truncating the action at the leading order in the lattice spacing a. Let us also note that an analogous 
model with 9 (instead of 3) Hermitian matrices has been studied by Monte Carlo simulation [35, 51] from 
different motivations. In that case the phase transition occurs at T ~ 0.9, which is slightly lower than the 
present model. 

n Let us note that the results in the fuzzy sphere phase are also consistent with Eguchi-Kawai's statement. 
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the volume independence of single-trace operators in D-dimensional U(oo) gauge theory 
provided that the U(l) z? symmetry is not spontaneously broken [53]. 

We have performed a similar analysis at a = 1.8, which is barely below the boundary 
of the fuzzy sphere phase. (See figure 2.) The Hagedorn temperature turned out to be 
Th ~ 1.1 as well. In the Yang- Mills phase, the Chern-Simons term M takes small values as 
one can see from figure 1, and the observables have little dependence on a. This property 
is found also in the totally reduced model studied in ref. [48] . Note also that the Hagedorn 
temperature Th is an 0(1) quantity, which means that Th = iV _2 / 3 TH vanishes in the 
N — > oo limit. In other words, if we drew the critical line corresponding to the Hagedorn 
transition in figure 2, it would be pushed towards the T = line in the large- iV limit. This 
is simply a reflection of the fact that, in the fuzzy sphere phase, we have to consider super 
high temperature to see non-trivial temperature dependence. 

6. Fuzzy-sphere/ Yang-Mills transition at high temperature 

In general, field theories at high temperature are effectively described by bosonic field 
theories in one dimension less. This phenomenon provides a useful approach to QCD at 
high temperature. (See, for instance, refs. [54,55] and references therein.) In the present 
model, 12 we do not have the subtlety related to infrared divergences unlike in ordinary field 
theories, since there is no infinitely extended spatial directions from the outset. 

The dimensionally reduced model is obtained from the original action (2.1) by sup- 
pressing the t dependence of the Id fields as 

S DR = ^{-^([AX,]) 2 - ^([X^X^ 2 + hae ljk X t X 1 X k } 

= N^-^A^A^Y + h 1 e ijk A i A j A k } , (6.1) 
where the Greek indices /x, v run over 1, • • • ,4 and we have defined 

Ai = T" 1 / 4 ^ (i = 1, 2, 3) , A A = T-^A , (6.2) 
7 = T^a . (6.3) 

The observables studied in the previous sections can be obtained at high temperature as 

{R 2 )c , T i/i./l tl{Ai f\ , ( 6 .4) 



N I DR>7 

(M) ~ T 3 / 4 • ( ^1 e ijk triAAjAk)) , (6.5) 

' ' 1 DR, 7 

/- ^ -r-(U:{[A,. I,]) ) , (6.6) 

DR, 7 



<|P|>^1- VWltr^) 2 / . (6.7) 



2^ \N / 



DR, 7 



If we fix T rather than T in the large- iV limit, the Polyakov line vanishes identically, and all the observables 
have no dependence on T. On the other hand, if we fix T in the large- N limit, the Polyakov line vanishes 
only at T = 0, and all the observables have non-trivial dependence on T. 

12 The high temperature limit in matrix quantum mechanics is also discussed refs. [35,56,57]. 
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The symbol ( • }dr,7 represents the expectation value with respect to the dimensionally 
reduced model (6.1), where 7 is related to a through (6.3). 

In the a = case, the corresponding dimensionally reduced model (6.1) is studied in 
detail at large N [58]. For instance, we have 



-/i, r ([A„^]) ; 



/ DR,0 



2.162(5) , 
1 



DR,0 



4 1- 



N 2 



(6.8) 
(6.9) 



Taking into account that the Greek indices run from 1 to 4 in contrast to the Roman 
indices, which run from 1 to 3, we obtain the asymptotic behavior of the original model 
with a = at high T as 



lim (R 2 ) ~ -CVT , 

N^oo 4 



(6.10) 
(6.11) 



Figures 4 and 5 show that our Monte Carlo results approach these results at high T. (Small 
deviations can be nicely reproduced by the next-leading order calculation [57].) 

In the fuzzy sphere phase, we can confirm the dimensional reduction analytically by 
using the all order calculation in perturbation theory. By taking the (3 — ► limit in the 
results for the full model, we obtain the all order results for the dimensionally reduced 
model, which can be obtained similarly to ref. [46]. In contrast to the situation in the 
Yang-Mills phase, T instead of T has to be large (in the large-iV limit) in order for the 
dimensional reduction to take place. 
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Figure 6: (Left) The observable (-^ tr (Aj) 2 ) DR in the dimensionally reduced model is plotted 
against 7 for TV = 12, 16, 18. The open and closed symbols represent the results for the single fuzzy 
sphere start and the zero start, respectively. The solid lines represent the all order results. (Right) 
The upper and lower critical points represented by closed and open circles, respectively, are plotted 
against N in the log-log scale. The straight lines represent the fits to 7^ = c\ and 7$ = c 2 TV _1 / 2 , 
where a = 0.9765 and c 2 = 2.5160. 
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Using the dimensionally reduced model (6.1), let us investigate the phase transition 
between the fuzzy sphere phase and the Yang-Mills phase in the high temperature limit. 
This clarifies, in particular, the first order nature of the phase transition, and it also enables 
us to make explicit the connection to the known results in a totally reduced model [48]. 
We perform Monte Carlo simulation 13 of the dimensionally reduced model (6.1) using, as 
the initial configuration, either of the two configurations given by 

_ j 7 (the single fuzzy sphere start) , ^ ^ 

1 (the zero start) , 

and ^4 = for both cases. In figure 6 (Left) we plot the observable appearing on the right 
hand side of eq. (6.4) against 7 for N = 12, 16, 18. For comparison we also plot the all 
order results for the dimensionally reduced model obtained from the perturbation theory 
around the single fuzzy sphere in the large-iV limit. The Monte Carlo results depend on 
the initial configuration in the intermediate region of 7, and we observe discontinuities at 

f 7c? ~ -7= for the single fuzzy sphere start , 
[ 7cr ~ 0.98 for the zero start , 

which we call the lower/upper critical points, respectively. (See figure 6 (Right) for a plot 
showing the large- N behaviors.) This clearly demonstrates that the phase transition is of 
first order. 

In a similar model [48], which can be obtained by simply omitting A4 from (6.1), the 
critical points are obtained as 7$ ~ and 7^ ~ 0.66, respectively. We find that the 
inclusion of the fourth matrix A4 changes the numerical coefficients, but not the powers of 
N, in the large- N behavior of the critical points. 

Using the relation (6.3), we obtain the critical points in terms of the parameters of the 
full model as 



a. 



T l/4 ? 



«= { f \ VN , (6-14) 

In this terminology, the critical point a c shown in figure 2 is actually the lower critical 
point. Note that the factor -7= in (6.14) is absorbed by the rescaling (3.12) of a and T, 

and the result agrees with the high T behavior (3.16) obtained from the effective action. 



7. Summary and discussions 

We have studied thermodynamical properties of a fuzzy sphere in a BFSS-type matrix 
model including the Chern-Simons term. We have established the phase diagram in the 
(a, T)-plane, and obtained, in particular, the phase boundary between the fuzzy sphere 
phase and the Yang-Mills phase as shown in figure 2. 

In the fuzzy sphere phase, we are able to obtain all order results for various observ- 
ables exploiting the one- loop saturation of the effective action in the large- iV limit. This 

13 We have used the same algorithm as in ref. [48]. 



- 15 - 



technique was previously applied to various fuzzy manifolds in totally reduced models. We 
consider it interesting that it can be generalized to a finite temperature setup in a straight- 
forward manner. Following refs. [49, 50] thermodynamical properties of four-dimensional 
fuzzy manifolds such as fuzzy CP 2 and fuzzy S 2 x S 2 can be studied in a similar way. 

One of the interesting aspects of our results is the scaling of parameters in the large- iV 
limit. In the fuzzy sphere phase, if one fixes the original parameters a and T in the large- ./V 
limit, one simply obtains trivial results corresponding to the classical fuzzy sphere at zero 
temperature. In order to keep non-trivial quantum corrections and thermal effects, one has 
to hold a and T fixed in the large- N limit. 

In that limit, we find that the Polyakov line vanishes smoothly as T approaches 0. 
This implies that the fuzzy sphere phase is not further divided into the confined phase and 
the deconfined phase, unlike the Yang-Mills phase. If we take the large- N limit at fixed T, 
we are always in the confined phase. If we take the large- N limit at fixed nonzero T, we are 
always in the deconfined phase. In ref. [23] it is stated that the Hagedorn temperature for 
the fuzzy sphere is Th = oo in an analogous model. We consider that our results provide 
a more precise formulation of that statement. 

As an outlook, we note that fuzzy manifolds [59-61] are also studied intensively in the 
IIB matrix model [47] in order to investigate the dynamical generation of 4d space-time. 
The same issue has been addressed by various approaches [62-69], and in ref. [68] the first 
evidence for such a phenomenon is obtained by the Gaussian expansion method. Based 
on the Eguchi-Kawai equivalence [53], two of the authors (N.K. and J.N.) conjectured [52] 
that a similar phenomenon should occur in the BFSS matrix model [12]. We therefore 
consider that studying the effective action for fuzzy manifolds in the BFSS matrix model 
would be an interesting future direction. In that case, the effective action is expected to 
be saturated at two loop similarly to the situation in the IIB matrix model [45]. 

From the view point of the gauge/gravity correspondence, the fuzzy sphere solutions 
in the pp-wave matrix model can be interpreted as giant gravitons. It would be interesting 
to look for phenomena in the dual gravity theory corresponding to the ones discussed in 
this paper. 
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